Fidelity matters: the birth of entanglement in the mixing of Gaussian states 
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We address the interaction of two Gaussian states tlirougli bilinear exchange Hamiltonians and 
analyze the correlations exhibited by the resulting bipartite systems. We demonstrate that entangle- 
ment arises if and only if the fidelity between the two input Gaussian states falls under a threshold 
value depending only on their purities, first moments and on the strength of the coupling. Our result 
clarifies the role of quantum fluctuations (squeezing) as a prerequisite for entanglement generation 
and provides a tool to optimize the generation of entanglement in linear systems of interest for 
quantum technology. 
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Gaussian states (GS), that is quantum states with 
Gaussian Wigner functions, play a leading role in contin- 
uous variable quantum technology [l[ for their extremal 
properties and because they may be generated with 
current technology, in particular in the quantum optics 
context 0-IH- As a consequence, much attention have 
been dedicated to the characterization of Gaussian en- 
tanglement [6-14|. Among the possible mechanisms to 
generate Gaussian entanglement, the one consisting in 
mixing squeezed states |15l42lj is of special interest in 
view of its feasibility, which indeed had been crucial to 
achieve continuous variable teleportation p^ . The en- 
tangling power of bilinear interactions has been widely 
analyzed, either to optimize the generation of entangle- 
ment 2^ 24 1 or to find relations between their entangle- 
ment and purities [i^ or teleportation fidelity [26. 27|. 

In this Letter we address bilinear, energy conserving, 
i.e., exchange, interactions described by Hamiltonians of 
the form Hj = g{a^h + ahi^), where a and h are bosonic 
annihilation operators, [a,a^] = 1 and [b,h^\ = 1, and 
g the coupling constant. Hamiltonians of this kind are 
suitable to describe very different kinds of quantum sys- 
tems, such as, e.g., two light- mode in a beam splitter or 
a frequency converter, collective modes in a gases of cold 
atoms [i^l , atom- light nondemolition measurements [i^ , 
optomechanical oscillators 27 , [s^ , nanomechanical os- 
cillators 31 1, and superconducting resonators [32|, all of 
which are of interest for the quantum technology. Our 
analysis can be applied to all these systems and lead to 
very general results about the resources needed for Gaus- 
sian entanglement generation. 

The bilinear Hamiltonians Hj generally describe the 
action of simple passive interactions and, in view of this 
simplicity, their fundamental quantum properties arc of- 
ten overlooked. Actually, the exchange amplitudes for 
the quanta of one of the systems strongly depend on the 
statistics of the quanta of the other one and on the par- 
ticle indistinguishability. This mechanism gives rise to 
interference and, thus, to the birth of correlations in the 



output bipartite system. A question arises about the 
nature of these correlations, depending on the parame- 
ters of the input signals and coupling constant. In this 
Letter, motivated by recent results on the dynamics of 



bipartite GS through bilinear interactions [33j, |3j| and 
by their experimental demonstration [ssf . we investigate 
the relation between the properties of two input GS and 
the correlations exhibited by the output state. Our main 
result is that entanglement arises if and only if the fi- 
delity between the two input states falls under a thresh- 
old value depending only on their purities, first-moment 
values and on the strength of the coupling. Our analysis 
provides a direct link between the mismatch in the quan- 
tum properties of the input signals and the creation of 
entanglement, thus providing a better understanding of 
the process leading to the generation of nonclassical cor- 
relations. In fact, if, from the one hand, it is well known 
that squeezing is a necessary resource to create entangle- 
ment [16|, |l7|, [24 1 , from the other hand in this Letter we 
show what is the actual role played by the squeezing, that 
is making the two input GS different enough to entangle 
the output system. 

The most general single-mode Gaussian state can be 
written a.s g = g{a,(, N) = D{a)S{^)vth{N)S^ {()D^ (a), 
where S{r) = cxp[i(^a^^ — ^*a^)] and D{a) = exp[aa'l' — 
a* a)] are the squeezing operator and the displace- 
ment operator, respectively, and i'th{N) = (N)'^ °/(l -|- 
^-jata-i-i thermal equilibrium state with N average 
number of quanta, a being the annihilation operator. Up 
to introducing the vector of operator = (i?i,i?2) = 
{q,p), where g = (a + a^)/\^ and p ~ {a — a^)/(i^/2) are 
the so-called quadrature operators, we can fully char- 

T 

acterize g by means of the first moment vector X = 
(R^) = y2(5Re[a],9m[a]), with (A) = Tr[Ag], and 
of the 2x2 covariance matrix (CM) cr, with [(t];,^. = 
^{RhRk + RkRh) - {Rh){Rk). k = l,2, which explicitly 
reads: [(T]kk = (2^)""^ [cosh(2r) — (— l)*-' cos(-0) sinh(2r)] 
for fc = 1,2 and [(t]i2 = M21 = -(2^)-i sin(V') sinh(2r). 
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rity of the state fj, = Tr[£<^] = (1 + 2A'')~^. Since we are 
interested in the dynamics of the correlations, which are 
not affected by the first moment, we start addressing GS 
with zero first moments (a = 0). The general case will 
be considered later on in this Letter. 

When two uncorrelated, single-mode GS Qk with CMs 
(Tj,, k = 1,2, interacts through the bilinear Hamilto- 
nian Hj, the 4x4 CM S of the evolved bipartite state 
£'12 = Ug{t)gi (K) Q2Ul{t), Ug{t) = ex^{-iHjt} being 
the evolution operator, can be written in the following 
block- matrix form [H : 



Si = T(Ti -f (1 - t)(T2, 
^2=T(T2 + {1-t)(Ti, (1) 
Sl2 = t(1 - t){(T2 - CTi), 



T = cos^ (gt) being an effective coupling parameter, and 
where the presence of a nonzero covariance term S12 sug- 
gests the emergence of quantum or classical correlations 
between the two systems. Since S12 depends on the dif- 
ference between the input states CMs, a question natu- 
rally arises about the relation between the "similarity" 
of the inputs and the birth of (nonlocal) correlations. 
In this Letter we answer this question and demonstrate 
that entanglement arises if and only if the fidelity be- 
tween the two input GS falls under a threshold value, 
which depends only on their purities, the value of the 
first moments, and the coupling r. 

Let us now consider the pair of uncorrelated, single- 
mode GS Qk = Qi^k,Nk), fc = 1,2 and assume, without 
loss of generality, ^1 = ri and ^2 = r2e'^^ , with rk,tp G R. 
After the interaction, we found that the presence of en- 
tanglement at the output is governed by the sole fidelity 
Q2) = [Tr(-y^7g7e'2A/gi)]^ between the inputs. Our 
results may be summarized by the following: 

Theorem 1 The state gi2 ~ Ug{t) gi (E) g2 U^{t), result- 
ing from the mixing of two GS with zero first moments, 
gi{ri,Ni) and g2{r2e'^^ , N2) , is entangled if and only 
if the fidelity F{gi,g2) between the inputs falls below a 
threshold value i^c(/^i, IJ-2tT), which depends only on their 
purities fik — Tr[gf.] = (1 -|- 2Nk)~^ , k = 1, 2, and on the 
effective coupling parameter r = cos^ (gt) . 

Proof: In order to prove the theorem, we recall that a 
bipartite Gaussian state gi2 is entangled if and only if the 
minimum symplectic eigenvalue A of CM associated with 
the partially transposed state is A < 1/2 Q. Moreover, 
without loss of generality, we can address the scenario in 
which rk and Nk, k = 1,2, are fixed and we let tp vary in 
the interval [0, 27r]. 

First of all, we prove that A < 1/2 F{gi,g2) < 
fc(Mij /^2; t) (necessary condition). As we will see, this 
will allow us to find the analytic expression of the thresh- 
old Fe(/ii, /i2; r), which will be used to prove the suffi- 
cient condition, i.e., F{gi,g2) < Fc{^i, H2',t) ^ A < 1/2. 
Fig. [H shows the typical behavior of A and of the fidelity 




FIG. 1: (Color online) Plot of the fidelity F(gi, Q2) (red, solid 
line) between the two input states and of the minimum sym- 
plectic eigenvalue A as a function of iji for r = 0.5 (blue, 
dashed line) and r = 0.8 (purple, dotted line). The other 
involved parameters are ^1 = 0.5, A^i = 0.2, ^2 ~ 0.7e** and 
= 0.3. The colored regions denote the ranges of leading 
to an entangled state for the given r, while the horizontal 
dot-dashed lines refer to the corresponding thresholds Fc- 



F as. & function of the squeezing phase ip for fixed r^ and 
k = 1,2 (here we do not report their analytic expres- 
sions since they are quite cumbersome). As one can see, 
both A and F are monotone, decreasing (increasing) func- 
tions of ip in the interval [0, tt) ([tt, 27r], respectively) and 
have a minimum at tJj ~ n, whose actual value depends 
on both Tfc and Nk but not on r. In our case, one finds 
that, for fixed r^ and N^, k = 1,2, if min^ A < 1/2, then 
there exists a threshold value ip^ = V'o(^ij Mij M2j ''"): 



arccos 



J' cosh(2ri)cosh(2r2) - / (^1,^2, 
\ sinh(2ri) sinh(2r2) 



where we introduced: 



/(Ml,Ai2,T) 



l + /i?Mi-(A^?+Mi)(l-2r)2 



*A*iM2t(1 - r) 



and fik = Tr[£»|] = (1 + 2Nk)''^, k = 1,2, are the pu- 
rities of the inputs, such that ii ip £ {ipig, 27r — -00) then 
A < 1/2, i.e., gi2 is entangled. Since the fidelity be- 
tween the two GS gk, characterized by the CMs cr^, 
k = 1,2 (and zero first moments), is given by [36 1 

Figi, g2) = (VA + 5- , where A = det[cri + (T2] 

and 5 ~ 4 J|^.^j^(det [crj.] — |), the threshold value F^ = 
Fo(/ii, /X2; t) of the fidelity is thus obtained by setting 
ip — tpc and explicitly reads: 



F, = 



4mi/^2V'^(1 - r) 



v/.g- + 4r(l - T)g+ - v/4t(1 - r) 



r)g- 



(2) 



where g± = g±{ni,ii2) = nfc=i,2(-'^ ^ A*fc)- ^he threshold 
depends only on r and on the purities /Xfc of the input 
GS and is independent of the squeezing parameters r^, 
despite the fact tpe does. Finally, if A < 1/2, i.e., £112 is 
entangled, then F[gi, P2) < Fc{fii, ^2] t). This concludes 
the first part of the proof. 
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Now we focus on the sufficient condition, i.e., 
F{qi,Q2) < Foiiii, ii2]t) =^ \ < 1/2. Thanks to the 
first part of the theorem and since both F and A are 
continuous functions of tp, for fixed r^. and iV^., fe = 1, 2, 
which have a minimum in "0 = tt, it is enough to show 
that Fmin = Txiiujj, F(gi, g2) < Fo(/ii, t) ^ Anii„ = 
min^ A < 1/2. Wc have: 



F ■ — 

^ mm — 



^1 + /i2 + 2^i^2 cosh[2(ri + r2)] 
where g_ is the same as in Eq. ([2]), and: 



2 V^Hifi2 

with7 (/^^+^^)(l-2T)2+8/ii^2T(l-T)cosh[2(ri+r2)], 
respectively. The inequality Fmin < Fc{^i, fi2',T), where 
Fc(/ii, ^2; t) is given in Eq. ([2|), is satisfied if 7 > l+fiffi^^ 
which leads to Amin < 1/2, as one may verify after a 
straightforward calculation. Now, since A is a continuous 
function of ip, there exists a range of values centered at 
■0 = TT, where the minimum occurs, in which A < 1/2 and, 
thus, F{qi,Q2) < Fc{fii, fJ.2',T), because of the first part 
of the theorem (necessary condition) . This concludes the 
proof of the Theorem. □ 

As a matter of fact, the presence of nonzero first mo- 
ments does not affect the nonclassical correlations exhib- 
ited by a bipartite Gaussian state, which only depend on 
the CM p. Thus, we can state the following straightfor- 
ward: 

T 

Corollary 1 If = Tr[{qk,Pk) Qk] ^ 0, where qt = 
(afc -|- al.)/\/2 and p = (aj, — a^)/ (iV^) are the quadrature 
operators of the system k ~ 1,2, then the state Q12 = 
Ug (t) Qi ® Q2 U^g (t) is entangled if and only if: 



F{qi, Q2) < r(Xi,X2)Fo(^i,Ai2;T), 
where Fdpi, fJ-2',T) is still given in Eq. (0) and: 

r(Xi,X2) = CXp Uixf2('Ti +,T2)-lXi 



(3) 



(4) 



where X 
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(^1 



X, 



Proof: The proof follows from Theorem [T] by not- 
ing that the presence of nonzero first moments does 
not modify the evolution of the CM, whereas the 
expression of the fidelity becomes [s^ F{qi,Q2) = 

r(Xi, X2) A + (5 — \f5^ , where A and 5 have been 
defined above. □ 

Theorem [1] states that if the two Gaussian inputs are 
"too similar" the correlations induced by the interaction 
are local, i.e., may be mimicked by local operations per- 
formed on each of the systems. The extreme case cor- 
responds to mix a pair of identical GS: in this case the 



interaction produces no effectj_since the output state is 
identical to the input one [SS, 34|, i.e., a factorized state 
made of two copies of the same input states, and we have 
no correlations at all at the output. Notice that for pure 
(zero mean) states the threshold on fidelity reduces to 
i^c(l, 1; ■T') = 1 Vr, namely, any pair of not identical (zero 
mean) pure GS gives raise to entanglement at the out- 
put. On the contrary, two thermal states Vk = vt\i[Nk), 
k = 1,2, as inputs, i.e., the most classical GS, lead to 
F{vi,V2) > Fe{ni, ^j,2',t): this fact, thanks to the The- 
orem [1] shows that we need to squeeze one or both of 
the classical inputs in order to make the states different 
enough to give rise to entanglement. Notice, finally, that 
the thresholds in Eqs. ([2]) and ([3]) involve strict inequal- 
ities and when the fidelity between the inputs is exactly 
equal to the threshold the output state is separable. 

The threshold Fe{pii, fj.2] t) is symmetric under the ex- 
change fii <-> fi2 and if one of the two state is pure, i.e., 
if fj.k = 1 then Feir) = V^idh/ y/1 + fif, with /i 7^ fc, i.e., 
the threshold no longer depends on r. 

For what concerns Gaussian entanglement, i.e., the 
resource characterized by the violation of Simon's con- 
dition on CM 0, our results also apply to the case of 
non-Gaussian input signals, upon evaluating the fidelity 
between the GS with the same CMs of the non-Gaussian 
ones. In fact, violation of the Simon's condition is gov- 
erned only by the behavior of the CM independently on 
the Gaussian character of the inputs states. On the other 
hand, identical non-Gaussian states may give raise to en- 
tangled output, the mixing of two single-photon states in 
quantum optical systems being the paradigmatic exam- 
ple (STj . In other words, the entanglement raising from 
the mixing of two identical non-Gaussian states cannot 
be detected by Simon's condition on CM. 

Up to now we have considered the correlation proper- 
ties of the output states with respect to the fidelity be- 
tween the input ones. However, similar relations may 
be found for the fidelities F{gh,Qk), k,h = 1,2, be- 
tween the input and output states, respectively, where 
Qh = Trfc[pi2], with h ^ k, are the reduced density 
matrices of the output states taken separately. In this 
case we found that the output is entangled if and only if 
F{Qh,Qk) < Fe(£'/i,^fe) where all the thresholds j;(£i;i,efc) 
still depends only on ni, 112 and r (here we put as ar- 
guments the density matrices in order to avoid confusion 
with the previous thresholds). The analytic expressions 
of Fc{gh,Qk) are cumbersome and are not reported ex- 
plicitly, but we plot in Fig. [2] the input-output fidelities 
and the corresponding thresholds for a particular choice 
of the involved parameters. If we look at the interaction 
between the two systems as a quantum noisy channel for 
one of the two, namely, gk — >■ £{gk) = Tr/i[pi2], h ^ k, 
than the birth of the correlations between the outgoing 
systems corresponds to a reduction of the input-output 
fidelity: the correlations arise at the expense of the infor- 
mation contained in the input signals. In turn, this result 
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may be exploited for decoherence control and preserva- 
tion of entanglement using bath engineering [35| . 




FIG. 2: (Color online) Plot of the fidelities F{gh, Qk) for 
r = 0.8 and the same choice of the other involved parameters 
as in Fig. [T] The yellow region shows the interval of values of 
t/j leading to an entangled state. The right panel is a magni- 
fication of the green, boxed region of the left panel: the hor- 
izontal lines refer to the corresponding thresholds Fc{Qh, Qk)- 
See the text for details. 

In conclusion, we have analyzed the correlations ex- 
hibited by two initially uncorrelated GS which interact 
through a bilinear exchange Hamiltonian. We found that 
entanglement arises if and only if the fidelity between the 
two inputs falls under a threshold value depending only 
on their purities, the first moments, and on the coupling 
constant. Similar relations have been obtained for the 
input-output fidelities. Our theorems clarify the role of 
squeezing as a prerequisite to obtain entanglement out 
of bilinear interactions, and provide a tool to optimize 
the generation of entanglement by passive (energy con- 
serving) devices. Our results represent a progress for the 
fundamental understanding of nonclassical correlations 
in continuous variable systems and may found practical 
applications in quantum technology. Due to the recent 
advancement in the generation and manipulation of GS, 
we foresee experimental implementations in optomechan- 
ical and quantum optical systems. 
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